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Abstract 



We obtain a new phantom black plane solution in 4D of the Einstein- 
Maxwell theory coupled with a cosmo logical constant. We analyse their 
basic properties and obtain the extensive and intensive thermodynamic vari- 
ables, as well as the specific heat and the first law. Through the specific 
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heat and the so-called geometric methods, we analyse in detail their ther- 
modynamic properties, the extreme and phase transition limits, as well as 
the local and global stabilities of the system. The normal case is shown with 
an extreme limit and the phantom one with a phase transition only for null 
mass. The systems present local and global stabilities for certain values of 
the entropy density with respect to the electric charge, for the canonical and 
grand canonical ensembles. 

Pacs numbers: 04.70.-s; 04.20. Jb; 04.70.Dy. 

1 Introduction 

It is well known that a black hole can radiate a black-body radiation 
when one takes into account the effects of classical gravitational field on 
quantized matter fields, i.e, a semi-classical analysis of the gravity [1]. So, 
we can make a study of the thermodynamic system of each new black hole 
solution. The most common method in the literature is the analysis made 
through the specific heat of the black hole [2] , which informs us if the system 
is thermodynamically interacting, if there exists any case in which the black 
hole is extreme or it passes across a second order phase transition. 

Recently, attention is attached to the methods for analysing the thermo- 
dynamic system through the geometry of the so-called thermodynamic space 
of the equilibrium states. The most common are the methods of Weinhold 
[3], Ruppeiner [4], geometrothermodynamics [5] and that of Liu-Lu-Luo-Shao 
[6] . These methods also notify if the system possesses thermodynamic inter- 
action and if it undergoes a second order phase transition, in addition to the 
properties about the stability. 

In this work, we desire to make a detailed analysis of the thermodynamic 
system of a well known class of solutions, with a particularly interesting 
symmetry, the planar. This class of solutions has been previously obtained 
for the case of planar and static symmetry in 4D, by Cai and Zhang [7]. 
This symmetry was then applied to traversable wormholes [8], and later, 
generalized to topological black holes in [9], and its various applications. We 
focus our attention to a class of solutions, called phantom [10], but now with 
a planar symmetry. 

Before beginning the analysis of this new class of phantom black holes, 
we will present briefly our interest in obtaining and studying such exotic 
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solutions. With the discovery of the acceleration of the universe, various ob- 
servational programs of studying the evolution of our universe were deployed, 
including the relationship of the magnitude-versus-redshift type supernovae 
la and the spectrum of the anisotropy of the cosmic microwave background. 
These programs promote an accelerated expansion of our universe, which 
should be dominated by an exotic fluid and should have a negative pressure. 
Moreover, these observations show that this fluid can be phantom, i.e, with 
the contribution of the energy density of dark energy [11]. 

As the interest in obtaining these classes has increased, we also found 
ourselves wanting to analyse a specific phantom model. We can mention 
here some recent results in the literature, such as the wormhole solutions and 
conformal continuation [ ], the black hole solutions of Einstein-Maxwell- 
Dilaton theory, [13], the higher-dimensional black holes by Gao and Zhang 
[14], and the higher-dimensional black branes by Grojean et al [15]. Analysis 
were also made in algebraic structures of this type of phantom system, as 
the case of the algebra generated by metrics depending on two temporal 
coordinates, with D > 5, which provides phantom fields in 4D, fulfilled by 
Hull [16], and Sigma models by Clement et al [17]. Here, we will obtain and 
study the thermodynamic properties of a solution arising from the coupling 
of Einstein-Hilbert action with a field of spin 1, which can be Maxwell or 
anti-Maxwell (phantom), and a cosmological constant, where the spacetime 
possesses planar symmetry. 

The paper is organized as follows. In Section 2, we present a new phan- 
tom black plane solution. The causal structure of the solutions are studied 
and the thermodynamic variables are obtained. The first law of thermody- 
namics is established and the specific heat is calculated. In Section 3, we 
minutely study the thermodynamics of normal and phantom solutions, using 
the analysis through the specific heat, subsection 3.1, and through the geo- 
metric methods of Weinhold, subsection 3.2, the geometrothermodynamics, 
subsection 3.3, and that of Liu-Lu-Luo-Shao, subsection 3.4. We finish the 
section with the study of local and global stabilities in subsection 3.5. The 
conclusion is presented in Section 4. 



3 



2 The field equations and the black holes 
solutions 



The action of the theory is given by: 

S = J dx A ^g~ [K + 2nF> a/ F ia> + 2A] , (2.1) 

where the first term is that of Einstein-Hilbert, the second is the coupling of 
(anti) Maxwell field F^ v = d^A^ — d u A^ with the gravitation, and the third 
is the cosmological constant. Making the functional variation of the action 
(2.1) with respect to the field and the inverse of the metric, g^, using 
R = — 4A, we get the following equations of motion 

V M [FH = 0, (2.2) 
R»u = 2 V {^g^F 2 - F^F V A - A 9fMU . (2.3) 

Let us write the static and plane symmetric line element as 

dS 2 = A{r)dt 2 - B{r)dr 2 - C{r){dx 2 + dy 2 ) , (2.4) 

with r = \z\. We will also assume that the Maxwell field is purely electric 
and only depends on r. With (2.4), one can integrate (2.2) and obtain 

with q a real integration constant. Substituting (2.5) into the equations of 
motion (2.3), we obtain the equations 

A" 1 (A'\ 2 A'B' A'C / q 2 A , , 

I-2(l) ~2AB + Hc = 2B ("^- A )' < 2 - 6 > 

A " lfA ' V *» ^-^-(%)=™L£-a), (2.7) 



A 2 \ A J 2AB C BC \C J \'C 2 

A'C C" B'C Sf Q 2 , 1 



2AC C 2BC yC 2 

where the "prime" denotes the derivative with respect to r. Choosing the 
coordinates such that 

A(r) = B-\t) , C(r) = a 2 r 2 (2.9) 
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with A = —3a 2 , the solution of the equations of motion (2.6)-(2.8) is given 
by 



dS 2 = A(r)dt 2 - A-\r)dr 2 - C(r)(dx 2 + dy 2 ) , 

F = -£r } dr A dt , A{r) = a 2 r 2 - f + V ^ , C{r) = a 2 r 2 



(2.10) 



where m is the mass and q the electric charge of the (phantom) black plane. 
This is the same solution as that of [ ], for r\ — 1, and phantom black plane 
solution for rj = —1, obtained for the first time here. 

We can rewrite the solution in terms of the densities of mass M and 
electric charge Q, as calculated in [7], yielding 



dS 2 = A{r)dt 2 - A-\r)dr 2 - C{r){dx 2 + dy 2 ) 



idr A dt , A(r) = a 2 r 2 - + r/^g- , C(r) = a 2 r 



„2»,2 



C(r 



(2.11) 



One can calculate the horizon of this solution, vanishing A(r), obtaining 

AttM Att 2 Q 2 



2 2 

a r 



a 2 r 



+ rj- 



a^r 2 



0. 



(2.12) 



This solution possesses two complex and two real roots. The real roots are 
given by 



r± 



2k ± 



8vrM 



- Ik 



2k 



(2.13) 



\ 



7TM 



a' 1 



+ 



'ttM 



a 1 



/4tt 2 Q 5 
V 3a 6 



+ 



\ 



7TM 



a' 




V 



3a 6 



(2.14) 



For the normal solution, rj — 1, one has < r_ < r + , and for rj = — 1, 
the corresponding is r__ < < r+, with r + > |r_|. We observe that in the 
phantom solution, r_ is in the negative part, but here something happens 
that we do not have in the spherical symmetry, because as r± = [zi^l, one 
gets Zi(-i-) = ±r + and z 2 (±) = ±r_. As r_ < 0, one gets zi(~) < z 2 (+) < 
< Z2{-) < Then, the singular plan z = r s = z s = is covered 
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by the plans z = = z 2 ( + ),z = z 2 (-) an d z = z%{+)- In the case 

of spherical symmetry, the internal horizon r_ could not be achieved, for a 
solution of non-degenerate horizon. Hence, here we have a drastic change in 
the causal structure of the phantom black plane solution, whose singularity 
is covered by two horizons in the positive part of z. This could not occur 
in the phantom solutions with spherical symmetry, where just one horizon 
covered the singularity. However, another unusual event happens, where we 
get two horizons but with the property of non existence of extreme case, i.e, 
these horizons can never be equal, when we consider only real values. 

Now, we are interested in the geometrical analysis representing semi- 
classical gravitational effects of the black hole solutions as mentioned before. 
By semi-classical we mean quantize the called matter fields, while the back- 
ground gravitational field is treated classically. Therefore, we will work with 
the semi-classic thermodynamics of black holes, studied first by Hawking [1], 
and further developed by many other authors [18]. 

There are several techniques to derive the Hawking temperature law. 
For example we can mention the Bogoliubov coefficients [19] and the energy- 
momentum tensor methods [2, 18], the euclidianization of the metric [20], the 
transmission and reflection coefficients [21, 22], the analysis of the anomaly 
term [23] , and the black hole superficial gravity [24] . Since all these methods 
have been proved to be equivalent [25], then we opt, without loss of generality, 
to calculate the Hawking temperature by the superficial gravity method. 

The surface gravity of a black plane is given by [7]: 



K 



2V-9oo9i 



(2.15) 



where r + is the event horizon radius, and the Hawking temperature is related 
with the surface gravity through the relationship [1, 24] 

T=-. (2.16) 

Z7T 



Then, for the black plane solution (2.11), we get the surface gravity (2.15) 

as 

and the Hawking temperature (2.16) in this case is : 



2 2vrM 4fQ^ 

K = a r + + ~IJ2 ~ v-7-r , (2-ir) 



T = — 
2tt 



2 2ttM An 2 Q 2 

a r + + ~^TY 4 3 
crri crri 



(2.18) 
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We define the entropy per unit of area of the black plane as two times the 
quarter of the horizon area 



1 cy^r 1 

S = 2x- 4 A = ^±, (2.19) 

where the factor 2 is due to the contribution of two planes z = ±r + . 

From (2.11), we can calculate the electric potential scalar at the horizon 

r 

A = F w (r')dr> =^-. (2.20) 
J r= r+ a 2 r + 

+oo 

Let us check the first law for the solution (2.11). Taking the differential 
of the mass, isolated from (2.12), of the electric charge and of the entropy 
(2.19), we get 

dM = ( ^— t - n^^] dr + + n^-dQ , dS = a 2 r + dr + , (2.21) 
\ An ot z r+ J ot z r + 

which satisfies the first law of thermodynamics 

dM = TdS + r]A dq . (2.22) 

Note that we introduced a compensating sign rj in (2.22) due to the contri- 
bution of the negative energy density, in the phantom case, the field of spin 
1, F^ u , which provides a work with an inverted sign in the first law. 

As we need to study the thermodynamic system through the geometric 
methods, we must first write the mass in terms of the entropy and the electric 
charge. We can do this by isolating the mass in (2.12) and then replace r + 
in terms of the entropy 5 , with the use of (2.19), which yields 

M(S, Q ) = a2S2+ £ Q2 , (2.23) 
nay lb 



where we have the conditions Q 2 < (3a 6 /47r 2 )(7rM/o; 4 ) 4//3 for rj = 1 [ j 
(real horizon in (2.13)) and Q 2 < (a 2 S 2 /7r 2 ) for rj = —1. We also write 
the temperature and the electric potential in terms of the entropy and the 



5 Wc take r + = \/2S/a, with the sign of a > 0. The negative sign of a can be considered 
taking r + = —\j2Sja. 
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electric charge. Taking (2.18) and (2.20), for r + in terms of the entropy, we 
get 



T(S,Q) = Ap=W (2 . 24) 

We can then calculate the specific heat by the expression 
r _(d_M\ _(d_M_\ l(&M\ _ 2S (3a 2 S 2 - r]n 2 Q 2 ) 



dT ) q \dSj q / \dS 2 J q 3 (a 2 S 2 + r]7T 2 Q 

We now have in hand the basic requirements to begin our analysis of the 
thermodynamic system of these solutions. In the next section we will study 
the specific heat (2.25) and through the four geometric methods, the ther- 
modynamic properties of these planar solutions. 



3 Thermodynamics of black plane 

In this section we will study in detail the thermodynamic properties of 
the planar solutions (2.11), both for normal and phantom cases. Through 
the specific heat and the curvature scalar of the thermodynamic spaces of 
the equilibrium states, we will examine whether there is an extreme case 
(only by the usual method), phase transition and finally, the local and global 
stabilities of the thermodynamic system. 



3.1 Analysis of specific heat 

Historically, the study of specific heat for revealing the thermodynamic 
properties was the first to be used [2] and has been called of usual method. 

Here, we have the expression of the specific heat (2.25), which, equating 
to zero, reveals the value of the entropy for which the solution is extreme, 
i.e, for S = S e = rtQy/rj/ay/S, which is real only for rj — 1. Therefore, there 
does not exist an extreme case for the phantom solution with rj = —1, as we 
had seen in its causal structure. 

Similarly, we can find the value of the entropy for which the system un- 
dergoes a second order phase transition, i.e, when the specific heat diverges. 
In this case the specific heat (2.25) diverges for S = S t = —iy/rjuQ/a, which 
shows that the normal case r] — 1 has no phase transition, while the phantom 
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case possesses a phase transition in S = S t . Note that this case is the specific 
value where the mass (2.23) vanishes. 

We will take the results of the study of specific heat as the basis for 
comparing with a geometric analysis of the thermodynamic system, through 
the four most popular methods in the literature. All these methods have 
in common the definition of a metric for the thermodynamic space of the 
equilibrium states, where the calculation of the curvature scalar of this metric 
reveals the existence or not of thermodynamic interaction, phase transition 
points, among other thermodynamic properties. Let us calculate this object 
with the aid of a mathematical software. 

In the next subsection we will analyse the thermodynamic system through 
the method of Weinhold. 

3.2 The Weinhold method 

Historically, Weinhold was one of the first to formulate a geometric de- 
scription applicable to a thermodynamic system. The method of Weinhold 
[3], as it is known, aims to define a metric for the thermodynamic space of 
the equilibrium states, through the mass (2.23) as thermodynamic potential. 
The metric constructed in this way provides a curvature scalar Rw, which, 
for this method can be interpreted as a function of extensive variables that 
shows the points of phase transition, when there exists, where the thermo- 
dynamic system goes by. Then, we define the metric of Weinhold as being 

f) 2 M f) 2 M f) 2 M 



3(a 2 S 2 + r]n 2 Q 2 ) ]o2 2 V ttQ ^ ] ^ V V2 
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dS' - -=^-dSdQ + ^dQ' . (3.26) 

Here we see that the curvature scalar Ry/ of this metric is identically zero, 
which prevents us of doing an analysis of the phase transition of the thermo- 
dynamic system. This result does not agree with the study of the specific 
heat. In the next subsection we will study the thermodynamics through the 
method of geometrothermodynamics. 

3.3 The Geometrothermodynamics method 

The Geometrothermodynamics (GTD) [5] makes use of differential geom- 
etry as a tool to represent the thermodynamics of physical systems. Let 
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us consider the (2n + l)-dimensiorial space T, whose coordinates are rep- 
resented by the thermodynamic potential $, the extensive variable E a and 
the intensive variables I a , where a — 1, ...,n. If the space T has a non de- 
generate metric Gab(Z c ), where Z c = {&,E a ,I a }, and the so called Gibbs 
1-form 6 = 6?$ — 5 a bI a dE b , with 5 a b the delta Kronecker; then, the structure 
(T, 0, G) is said to be a contact riemannian manifold if G A (c/G) n ^ is 
satisfied [ 26]. The space T is known as the thermodynamic phase space. We 
can define a n-dimensional subspace Ec T, with extensive coordinates E a , 
by the map ip : E -> T, with $ = $(£ a ), such that y>*(G) = 0. We call the 
space E the thermodynamic space of the equilibrium states. 

We can then define the metric of the thermodynamic space of the equi- 
librium states E, through the derivation of the thermodynamic potential and 
its extensive variables as [27] 

<•) = ( EC w) dEadEb > ( 3 - 27 ) 

which, by definition, is invariant under Legendre transformations. Through 
the metric (3.27), we can calculate the curvature scalar of the space E, which 
informs if the system passes by a phase transition, when the scalar diverges 
for some value of extensive coordinates. If the scalar is not zero, the system 
possesses thermodynamic interaction, i.e, the Hawking temperature is non 
null. 

Here, we will do the calculation of the metric of E, using the mass (2.23) 
as the thermodynamic potential, which provides 

dl o(M) ^253 dS + dQ ■ (3 - 28) 

The curvature scalar of this metric is given by 

_ 8aWS* {-WS* + 7VQ 2 ) 

9 (a 2 S 2 + r/7r2Q2)4 • ^ ZJ ) 

We get the value for which the scalar (3.29) diverges, which is given by 
St = —i^nQ/a, in agreement with the value obtained through the specific 
heat (2.25). This result is consistent with the specific heat, where we have 
found that the normal case has no phase transition and in the phantom case 
has one point of second order phase transition in S — St- 

In the next subsection we will see the analysis made by the geometric 
method of Liu-Lu-Luo-Shao. 
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3.4 The Liu-Lu-Luo-Shao method 



The geometric method of the analysis of the more recent thermodynamic 
system is that of Liu-Lu-Luo-Shao [ ], which defines a metric in the ther- 
modynamic space of the equilibrium states, based on the Hessian matrix of 
several free energy, the Helmholtz's one in our case, and which can be written 
as follows 

dl LLLS{F) = -dTdb + r]dA dq = -^db + ~ ~dg~ J q + T] ~dq~ q 

- S (» !S! + "'« , W t ,* i «'. (3.30) 



The curvature scalar of this metric is given by 

V2a 3 nS 5 / 2 
3 {a 2 S 2 + r]7r 2 Q 



Rllls = - • (3-31) 



Then, the analysis by this method shows that the normal case does not 
possess phase transition and the phantom case possess a transition phase at 
S = S t = —i^TrQ/a, which is in agreement with the specific heat. 

In the next subsection we will study the local and global stabilities of the 
black plane solutions. 



3.5 The local and global stability 

Let us now study the local and global stabilities of these solutions. 
Through the specific heat (2.25) and the temperature (2.24), one can see that 
in the normal case, r\ — 1, the system is locally stable for 3a 2 S 2 > tt 2 Q 2 , with 
C q ,T > 0, and unstable for the other values. In the phantom case, r] — — 1, 
the system presents a local stability for a 2 S 2 > 7r 2 Q 2 , with C q , T, M > 0. 
Defining the Gibbs's potential 

G = M-TS-,,A Q=-( a2 i 2 + % Q2 ), (3.32) 

we get that in the normal case, in the grand canonical ensemble, the system 
is globally stable for any values of S and Q, with G < 0, VS 1 , Q. But in 
the phantom case, the system is globally stable only if a 2 S 2 > tt 2 Q 2 , which 
agrees with the local stability of the specific heat. 
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In the canonical ensemble, we can define the Helmholtz free energy as 



which yields a globally stable system (F < 0), for the normal case, when 
a 2 S 2 > 3tc 2 Q 2 , and for the phantom case F < 0,WS,Q. 

4 Conclusion 

We obtained a new phantom black plane solutions in (2.11). We anal- 
ysed their basic geometric properties, obtaining the thermodynamic vari- 
ables, temperature (2.18), entropy density (2.19) and the electric potential 
(2.20). We established the first law of thermodynamics in (2.22) and calcu- 
lated the specific heat (2.25). 

We analysed the thermodynamic system through the study of the spe- 
cific heat and the geometric methods called Weinhold, the geometrothermo- 
dynamics and that of Liu-Lu-Luo-Shao. In the Weinhold's case, the space 
metric is not invariant under Legendre transformations, and thus cannot rec- 
oncile a good thermodynamic analysis, therefore, in general, this method 
cannot agree with that of specific heat. By the use of the geometrothermo- 
dynamics and the method of Liu-Luo-Shao, we obtain the same results as 
in the case of specific heat, which shows that these two geometric methods 
agree with the usual one. 

The summarized results are that the normal case possesses an extreme 
limit for S = S e = irQ^/a\/3, and the phantom case presents a phase 
transition point in S = St = —i^nQ/a, which represents a solution with 
mass (2.23) identically null. The interpretation of massless solutions has been 
presented in [21], but without any conclusion about its thermodynamics. 

The normal case presents locally stable thermodynamic system, for 3a 2 S 2 > 
n 2 Q 2 , and globally stable, in grand canonical ensemble, when G < 0,WS,Q, 
and in canonical ensemble for a 2 S 2 > 37r 2 Q 2 . On the other hand, the phan- 
tom case is locally stable when a 2 S 2 > tt 2 Q 2 , and globally stable, in grand 
canonical ensemble, when a 2 S 2 > tc 2 Q 2 , and in canonical ensemble, when 
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(3.33) 



F < 0,VS,Q. 
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